Breakup of ^^F on ^°^Pb near the Coulomb barrier 
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Angular distributions of oxygen produced in the breakup of ^'^F incident on a ^"^Pb target have 
been measured around the grazing angle at beam energies of 98 and 120 MeV. The data are domi- 
nated by the proton stripping mechanism and are well reproduced by dynamical calculations. The 
measured breakup cross section is approximately a factor of 3 less than that of fusion at 98 MeV. 
The influence of breakup on fusion is discussed. 
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I. INTRODUCTION 

The study of nuclear reactions near the Coulomb bar- 
rier involving loosely bound nuclei has received consider- 
able attention in recent years. This is primarily driven by 
the advent of radioactive ion beams M. It is frequently 
observed in stable beam experiments that the subbarrier 
fusion cross sections are enhanced over one-dimensional 
barrier penetration model predictions. The enhance- 
ment can be described by channel couplings where the 
interplay of the intrinsic degrees of freedom and reaction 
channels modify the single barrier to multiple barriers 
The barriers appearing at lower energies are responsible 
for the fusion enhancement. Breakup is a major reac- 
tion channel in the scattering of loosely bound nuclei; 
this removal of the incident flux would lead to fusion 
suppression ||[ |^. On the other hand, the coupling to 
the breakup channel can change the barrier distribution 
which could result in fusion enhancement Q . Experimen- 
tal efforts have been put forward to study the influence 
of breakup on subbarrier fusion 

The fusion excitation functions of the neutron skin nu- 
cleus ^Ee on 209Bi|, 0, | and ^^^V^ were measured 
and large subbarrier fusion enhancements were observed 
in both cases. The breakup cross sections of ^He on 
^''^Bi measured below the Coulomb barrier are orders 
of magnitude greater than fusion[^0[ 111]. Measurements 
with stable ^Be, which has a neutron binding energy of 
1.665 MeV, on ^"SpbQ and 209Bi|l|, |l| found that 
the fusion was not enhanced below the Coulomb barrier 



and was suppressed by about 30% above the barrier. In 
the ^Be+^^^Pb reaction, the incomplete fusion reaction 
Q,_l_208p|^ following the breakup of ^Be into n-l-a-l-Q: was 
observed. The suppression of fusion at energies above 
the barrier was attributed to the projectile breakup. 
The barrier distribution extracted from the fusion ex- 
citation functions is consistent with a single barrier for 
the ^Be-t-^^^Pb and ^Be-|-^°^Bi[|l5| even though very large 
breakup yields were observed below the barrier |l3]. The 
fusion of ^Be+^^^Bi was compared to that of ^^Be(a neu- 
tron halo nucleus with neutron binding energy of 0.504 
MeV)-|-^*'^Bi. At energies below the barrier the cross sec- 
tions were similar to that of ^Be+^^^Bi whereas at ener- 
gies above the barrier the cross sections are significantly 
larger than the predictions from a coupled-channels cal- 
culation which takes into account the large rms radius of 
^^Be. However, the precision of the data was not very 
good and further measurements are required p5[|. 

On the proton rich side, fusion of a proton drip line nu- 
cleus, ^^F, with ^°®Pb was measurcd|l^. The fusion ex- 
citation function is almost identical to that of ^^O+^^'^Pb 
and i9F-|-^°*Pb after correcting for the Coulomb barrier 
arising from the charge and size differences. There was 
no enhancement or perhaps a small suppression of fusion 
below the barrier. It is noted that the loosely bound 
proton can be polarized in the large Coulomb field of the 
target in such a way that the proton is shielded by the 
core and the breakup probability is reduced |l^. This 
paper reports the breakup of ^^F on ^osp]-, nieasured near 
the Coulomb barrier. 
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II. EXPERIMENT 

The experiment was carried out at the Hohfield Ra- 
dioactive Ion Beam FaciUty (HRIBF) where the Isotope 
Separator On-Line (ISOL) technique was employed for 
radioactive ion beam production. A 44 MeV deuteron 
beam from the Oak Ridge Isochronous Cyclotron (ORIC) 
was incident on a fibrous hafnium oxide target to produce 
short-lived by the i^O(d,n)i'^F reaction ||l9). The re- 
action products were extracted from a closely coupled 
kinetic ejection negative ion source pO|, mass analyzed, 
and accelerated by the 25 MV tandem electrostatic ac- 
celerator. The -'^^O isobar was removed from the accel- 
erated beams by inserting an 80 /ig/cm^ carbon foil at 
the exit of the tandem accelerator and selecting the 9+ 
ions with a 90 degree analyzing magnet. The target used 
was a self-supporting ^°*Pb foil with a nominal thick- 
ness of 1.8 mg/cm^. The reaction energies, 98 and 120 
MeV, were calculated for the beams at the middle of the 
target by taking into account the energy loss in the tar- 
get. The beam intensity was measured by detecting the 
secondary electrons generated during transmission of the 
beam through a 10 fig/cvo? carbon foil with a microchan- 
nel plate detector. The average intensity was 1.5 x 10^ 
and 8 x lO*' /s for the 98 and 120 MeV reactions, re- 

spectively. The highest intensity achieved was 10"^ ions/s 
for the 120 MeV beam. 

The reaction products were detected by a AE-E tele- 
scope composed of a 29 iim Si detector mounted in front 
of a 1000 fim double-sided silicon strip detector (DSSD). 
The area of the Si detector and the DSSD is 5 x 5 cm^ . 
The DSSD, which has 16 vertical and 16 horizontal strips, 
was placed near the grazing angle at 10.5 cm from the tar- 
get and symmetric with respect to the horizontal plane, 
i.e. half of the detector above and half of the detec- 
tor below the plane. At backward angles, the variation 
of scattering angle for pixels on a vertical strip is small. 
Events in pixels on the same vertical strip can be summed 
to increase statistics. At forward angles, pixels on the 
same vertical strip have to be divided into two groups, 8 
middle pixels and 8 outer pixels, in order to keep the an- 
gular spread similar to the backward angles (~ 2°). The 
uniformity of the Si detector was determined by measur- 
ing elastic scattering at forward angles. The position of 
the elastically scattered particles penetrating the detec- 
tor was obtained from the overlap of the horizontal and 
vertical strips of the DSSD. The energy loss of the elas- 
tically scattered particles in the 256 (16x16) pixels was 
compared to kinematics and stopping power calculations. 
Two 50 mm^ Si surface barrier detectors placed at 10° on 
either side of the beam were used to monitor the beam 
position and for normalization between runs. 

An E versus AE plot for the 120 MeV ^'^F-induced 
reaction is displayed in Fig. |^. It is obtained by summing 
events in the pixels of one vertical strip at 9iab — 65°. 
The energy loss, AE, was corrected for the nonuniformity 
of the Si detector. A group of oxygen events can be 
clearly identified and is well separated from the elastically 



scattered ^^F. 
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FIG. 1: Histogram of E versus A E for reaction products 
produced in 120 MeV "P+^^^Pb measured at Oiab = 65° by 
summing events in a vertical strip of the DSSD. 



III. RESULT AND DISCUSSION 

The angular distributions of oxygen produced in 
i7p_|_208p^3 coUisions at 98 and 120 MeV are shown in 
Fig. 1^ and Fig. ^, respectively. The angular distributions 
are bell shaped and have a peak near the grazing angle. 
Since the AE-E telescope is not able to resolve mass, 
calculations were performed to estimate contributions of 
reactions leading to oxygen isotopes other than ^^O. The 
charge exchange reaction (-'^'^F,^^0) has a Q-value of - 
0.11 MeV. Two-step Distorted Wave Born Approxima- 
tion (DWBA) calculations using the code FRESCO |2^ 
were performed to estimate the contribution of this reac- 
tion. Sequential single-nucleon transfer reactions, ^^F 

-> and i^F -> ^^P ^ "O, were calculated. In the 
calculations, transfer to excited states in the projectile- 
and target-like nuclei were included. States which have 
large spectroscopic factors measured in light ion trans- 
fer reactions or large cross sections calculated in one- 
step single-nucleon transfer reactions were selected for 
the two-step DWBA calculations. Table || presents the 
states included in the calculations. In these calculations, 
the spectroscopic factors were set to 1.0 to estimate the 
magnitude of the yields. The shape of the calculated 
(i^F,i'^0) angular distribution at Eiab = 120 MeV is sim- 
ilar to the measured angular distribution and has a peak 
at 0iab = 58°. However, the calculated peak cross sec- 
tion is 0.0028 mb/sr which is several orders of magnitude 
less than the measured value. Although reactions lead- 
ing to ^^O and ^"^Bi in the exit channel have positive Q- 
values, they cannot occur by simple single-step transfer 
processes. Therefore, the cross sections are expected to 

The re- 



be smaller than that for one-proton transfer] 22 



suits of DWBA calculations for one-step proton transfer 
208pb(i7p i6o)209Bj McV are shown by the dot- 
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ted curve in Fig. ||. One proton transfer to the six lowest 
single particle states in ^'^^Bi was calculated by the code 
Ptolemy [p3| with the spectroscopic factors set to 1.0. It 
can be seen that neither one-proton transfer nor charge 
exchange can account for the measured oxygen angular 
distribution. Since the direct charge exchange is orders of 
magnitude smaller than that of one-nucleon transfer p4[, 
its contribution to the data can be safely ignored. 
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FIG. 2: Angular distribution of oxygen produced from 98 
MeV ^^F-|-^°*Pb. The calculated stripping and diffraction 
breakup are shown by the dashed and dash-dotted curves, 
respectively. The solid curve is for the sum of the two. 
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FIG. 3: Angular distribution of oxygen produced from 120 
MeV ^''F+^°**Pb. The calculated stripping and diffraction 
breakup are shown by the dashed and dash-dotted curves, 
respectively. The solid curve is for the sum of the two. The 
results of one-step DWBA transfer calculations are shown by 
the dotted curve. 

The measured angular distributions are compared to 
results of dynamical calculations where the relative mo- 
tion of the proton and the ^^O core is described quantum 
mechanically by solving the time dependent Schrodinger 
equation for the two-body breakup in the Coulomb and 
nuclear fields from the target nucleus It has been 
shown that calculations of this kind are suitable for en- 
ergies near the Coulomb barrier p^. The time evolution 
of the projectile wave function was calculated to obtain 
the angular distribution of ^^O from the ^^F— >-'^^0-f p re- 
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action. The breakup angular distribution is obtained by 
multiplying the breakup probability calculated as a func- 
tion of impact parameter by a fit to the measured elastic 
scattering cross section at the corresponding Rutherford 
scattering angle. The measured angular distribution of 
elastic scattering at 120 MeV is shown in Fig. ^ Since 
the angular resolution is ~ 2° and the cross sections fall 
off exponentially at large angles, exponential functions 
fitted between two adjacent data points were used to cal- 
culate the weighted average of each data point at large 
angles. The solid curve shows an optical model fit to the 
elastic scattering, which is used in converting calculated 
breakup probabilities into an angular distribution. The 
resulting breakup cross section is shown by the dashed 
curve in Fig. ^ The dotted curve is the separate contri- 
bution from stripping, which dominates the breakup and 
is in fairly good agreement with the measurement shown 
by the open triangles. 

The calculated breakup cross section is compared to 
the measured angular distribution of oxygen fragments 
in Fig. ^ and Fig. |[ for 98 and 120 MeV, respectively. 
The dashed curve is the stripping and the dash-dotted 
curves is the diffraction dissociation, and their sum is 
shown by the solid curve. It can be seen that the mea- 
sured angular distribution is predominantly due to the 
stripping breakup reaction. The agreement between the 
data and calculations is very good for the 98 MeV mea- 
surement. For the 120 MeV reaction, the measured an- 
gular distribution is shifted slightly forward compared to 
the calculated distribution but the total cross sections 
are in good agreement (see Fig. ||). 

The angular distributions of oxygen were fit to a Gaus- 
sian function to obtain angle integrated breakup cross 
sections. Fig. ^ displays the fusion excitation function 
measured by Rehm et a/.|l6| and the breakup cross sec- 
tions measured in this work. The calculated diffraction 
and stripping breakup are shown by the dash-dotted and 
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FIG. 4: Angular distribution of elastic scattering (filled cir- 
cles) in 120 MeV ^^F+^°*Pb. The result of an optical model 
fit to the data is shown by the solid curve. The angular distri- 
bution of oxygen produced in the same reaction is presented 
for comparison (open triangles). The calculated stripping 
is shown by the dotted curve and the sum of stripping and 
diffraction breakup is shown by the dashed curve. 



dashed curves, respectively. As it was seen in the angu- 
lar distribution of oxygen fragments that the measured 
breakup is dominated by proton stripping, the angle inte- 
grated breakup cross sections are in good agreement with 
the calculated stripping cross sections. Near the barrier, 
the diffraction breakup is a factor of 3 less than stripping 
which, in turn, is about a factor of 3 less than fusion. 
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FIG. 5: Comparisons of the fusion excitation function (o) of 
i^p_j_208p|j measured by Rehm et a/, jlij and angle-integrated 
breakup cross sections measured in this work ( x ) . The dotted 
curve is the one-dimensional barrier penetration model predic- 
tion and the dashed and dash-dotted curves are for stripping 
and diffraction breakup, respectively, predicted by dynamical 
calculations. 

The result of a one-dimensional barrier penetra- 
tion model calculation for ^^F-|~^°^Pb using the code 
CCMOd[|^ is shown in Fig. || and Fig. || by the dot- 
ted curve. The barrier potential parameters, Vq — 
235.5 MeV, tq — 1.1 fm, and a — 0.65 fm, were taken 
from the analysis of fusion measurements of a neighbor- 
ing system, ^^0+^"®Pb[p8|, since the excitation function 



is almost identical to that of ^'^F+'^^^Ph after correcting 
for the Coulomb barrier. As can be seen, the calculation 
underpredicts the cross sections at subbarrier energies. 
Coupled-channels calculations were performed with the 
code CCMOD using procedures employed for analysis of 
the ^^0+^'^^Pb measurement in Ref. [g^. It is well estab- 
lished that inelastic excitations of the projectile and tar- 
get can contribute to subbarrier fusion enhancement. In 
many cases, coupling to inelastic excitation channels can 
account for the enhanced fusion rates. Calculations in- 
cluding the excitation of ^°^Pb to the lowest 2+, 3~, and 
5~ states were carried out. Furthermore, it was found 
in the analysis of the ^^0-|-^°*Pb data that the coupling 
of double-phonon excitations is essential for reproduc- 
ing the barrier distribution. The coupling of two-phonon 
states 3~ (8)3~ in the harmonic limit and all the resulting 
cross coupling terms, e.g. 3~ (8) 5~ were considered in the 
calculations. The result is shown by the dashed curve in 
Fig. ^. The calculation still underpredicts the measure- 
ment. The first excited state of ^''F is bound by 105 keV 
and can be excited from the ground state with a large 
B(E2) value, B(E2)i = 63.4 e^fm^g^^. The results of 
coupled-channels calculations including the excitation of 
^^F is shown by the solid curve in Fig. ^. It can be seen 
that the increase in the subbarrier cross sections is very 
small when this projectile excitation is included and the 
coupled-channels calculations still underpredict the sub- 
barrier cross sections. 
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FIG. 6: Fusion excitation function for "F-|-^°®Pb pre- 
dicted by a one-dimensional barrier penetration model (dot- 
ted curve), coupled-channels calculations taking into account 
target excitations (dashed curve), and projectile and target 
excitations (solid curve) as described in the text. 

In the fusion of ^^0+^°^Pb, the excitation function was 
reproduced by calculations coupling to the inelastic exci- 
tation channels only. Since all the Q-values for neutron 
transfer are negative, it is not necessary to consider trans- 
fer in the coupled-channels calculations. The neutron 
transfer Q-values in i'i'F-|-^''^Pb are positive for up to six- 
neutron pickup. In particular, the two- and four-neutron 
transfer have Q-values greater than 5 MeV. This can 
be compared to '"'Ca-l-^^Zr and *°Ca-|-^^Zr where very 
large subbarrier fusion enhancement was observed in the 
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latter [Q. Coupling to the inelastic excitations of pro- 
jectile and target only reproduces the ""^Ca+^^Zr mea- 
surement. There are still large discrepancies between the 
measured cross sections of ^°Ca-|-^^Zr and the coupled- 
channels calculations. The major difference in the two 
reactions is neutron transfer. The Q-values for multi- 
neutron transfer are negative in '*'^Ca-|-^'^Zr but positive 
in *'^Ca-f ®^Zr. Measurements of transfer near the barrier 
found large cross sections for the ^^Zr target The 
influence of transfer on fusion is demonstrated in a semi- 
classical model calculation!^ Q where the fusion 
cross sections as well as transfer for ''"Ca-l-^^Zr are re- 
produced. Fig. 1^ presents the results of coupled-channels 
calculations including the inelastic excitations discussed 
above and nucleon transfer. The transfer is treated ap- 
proximately in the code CCMOD, therefore only qualita- 
tive comparison can be made here. The transfer form 
factor is given by 



F{r) 



exp 



{r-Ri- R2) 



MeV, 



where Fq is the coupling constant, Ri,2 are the nuclear 
radii, and a=1.2 fm is the diffuseness parameter. The 
solid, dashed, dotted, and dash-dotted curves are for cou- 
pling constant Fq = 0.4, 1.0, 2.0, and 4.0 MeV, respec- 
tively. To simplify the calculation, three channels: one- 
proton stripping, one-neutron pickup, and two-neutron 
pickup, were included. The fusion excitation function can 
be reasonably reproduced with Fq set between 0.4 and 1 
MeV. However, it is noted that the quantity Fq can be as 
large as 3 or 4 MeV depending on the transferred angular 
momentum and the orbitals occupied by the transferred 
nucleons[27|, |3l|, Based on the calculations presented 
in Fig. the fusion excitation function can be repro- 
duced by including transfer of up to two nucleons in the 
calculations with the coupling constant Fq < 1 MeV. 
If channels of transferring more than two nucleons are 
included and Fq > 1 MeV is used, the calculation will 
overpredict the measured cross sections, i.e. the fusion is 
suppressed below the barrier. To better account for the 
influence of transfer on fusion in ^^F-t-^^^Pb, measure- 
ments of multi-nucleon transfer and more sophisticated 
model calculations such as in Ref. Q are required. Up to 
now, one of the reactions which has not been considered 
in the calculations is breakup. The simplified coupled- 
channels code used here cannot treat breakup rigorously. 
Nevertheless, it is conceivable that if fusion is suppressed, 
breakup can be responsible. 

It has been reported that the ^^F has a large rms ra- 
dius, < r >rms= 3.7 fm|3^. In the coupled-channels 
calculations, the nuclear radius is given by tqA^^^ where 
To = 1.1 fm is the radius parameter and A is the mass 
number. The effect of the large rms radius of ^^F was 
not accounted for in the previous calculations. To ex- 
plore these effects, the radius parameter of the projectile 
was adjusted in the calculations. Since the treatment 
of coupling to transfer degrees of freedom has large un- 
certainties introduced by the coupling constant Fq^ only 
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FIG. 7: Results of coupled-channels calculations including in- 
elastic excitations and transfer in ^"^F+^^^Ph. The results for 
transfer coupling constant Fo= 0.4, 1.0, 2.0, and 4.0 MeV are 
shown by the solid, dashed, dotted, and dash-dotted curves, 
respectively. 



inelastic excitations were included. In Fig. ^, the dashed, 
dash-dotted, and solid curves are for increasing the ra- 
dius of "F by 5 %, 10 %, and 20 %, respectively. It 
can be seen that the fusion excitation function can be 
well reproduced by increasing the radius of ^''F by 5 % 
whereas increasing the radius by 10 % results in a cal- 
culated cross section that exceeds the measurements. If 
transfer channels were included, the discrepancy would 
be larger. This further suggests that the fusion of ^^F 
and ^°^Pb may be suppressed. 
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FIG. 8: Results of coupled-channels calculations including 
only inelastic excitations (dotted curve) and the radius of the 
projectile increased by 5 % (dashed curve), 10 % (dash-dotted 
curve), and 20 % (solid curve). 

In the ^Be-f'^"*Pb reaction, complete fusion was found 
to be suppressed at energies above the barrier. How- 
ever, the sum of complete and incomplete fusion agrees 
with a coupled-channels calculation. The discrepancy 
between the measured complete fusion cross sections 
and the coupled-channels prediction is attributed to the 
breakup of ^Bc[^. The incomplete fusion arises from 
^Be breaking up into two a particles and a neutron, and 
subsequently an a particle fuses with the target. The 
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fusion measurements in Ref. [|6) are made by detecting 
fission fragments, and therefore probably determine total 
fusion-like cross section (complete + incomplete fusion) 
rather than the complete fusion cross section since the 
incomplete fusion reaction, ^^0+^"*Pb, produces fission 
events very similar to those of the complete fusion reac- 
tion. Consequently, it is not known whether the com- 
plete fusion of ^^F-|-^°*Pb is suppressed above the bar- 
rier. Measurements of ^^F on ^°^Pb at 10 MeV/nucleon 
showed that it is necessary to consider core absorption 
(^^O absorbed by ^°^Pb) in the dynamical calculation 
to reproduce the measured diffraction breakup yield [p 7 
The dynamical calculations presented in Fig. ^ and Fig, 
also include core absorption. It is expected that incom- 
plete fusion is present in ^^F-|-^°®Pb. A coincidence mea- 
surement of the breakup proton and the fission fragments 
are required to identify the incomplete fusion reaction. 

The predicted diffraction breakup seems too small to 
influence fusion significantly. However, the stripping 
breakup yield is about one third of fusion. The energy 
dependence of this reaction predicted by the dynamical 
calculation is presented by the dashed curve in Fig. ||. 
The calculated stripping cross section exceeds that for 
fusion below the barrier. The measured breakup of ^^^Li 
and ^Be in the vicinity of the barrier in ^'''Li-|-^'^^Pb|38 
and ^Be-|-^°^Bi[|l^, respectively, shows similar behav- 
ior. The analysis of elastic scattering of ^Li+^*'^Pb in 
Rcf. shows that the imaginary potential increases as 
the energy decreases below the barrier and the threshold 
anomaly[Q disappears. Because of this strong absorp- 
tion, the enhancement of fusion at low energies should 
be small and the breakup reaction is expected to be 
strong p^. The measured fusion yields for ^^F-l-^'^'^Pb 
were not enhanced and perhaps even slightly suppressed 
below the barrier. It is conceivable that strong absorp- 
tion exists resulting in large stripping breakup which re- 
moves ^^F from the fusion channel. It would be inter- 
esting to measure the elastic scattering and study the 
energy dependence of the interaction potentials. 

Large subbarrier fusion enhancement and trans- 
fer/breakup were observed in the ^He-t-^^^Bi reaction^ 
p^ . Analysis of elastic scattering indicated an absence 
of the threshold anomaly |pl|. In this case, the strong 
absorption may not enhance fusion much but may con- 
tribute mostly to transfer/breakup, as pointed out in 
Ref. . The neutron binding energy of ^He is fairly low. 
As suggested in Ref. 0, the large subbarrier fusion en- 
hancement may arise from neutron flow since the thresh- 
old barrier correlates with neutron binding energies pi]. 



In the -'^^F-|-^°^Pb reaction, the proton binding energy is 
very low but the proton flow must be strongly suppressed 
because of the Coulomb barrier. Therefore, the behav- 
ior of fusion below the barrier is different from that of 
6He+209Bi. 

IV. SUMMARY AND CONCLUSION 

The breakup of 98 and 120 MeV ^'^F on ^ospb was 
measured by detecting oxygen in singles. The angular 
distributions of oxygen are well reproduced by dynami- 
cal calculations and found to be dominated by stripping 
breakup. Near the barrier, the angle integrated stripping 
cross section is about 30% of that of fusion. It has been 
shown in the analysis of ^Li+-^"^Pb elastic scattering that 
the imaginary potential continues to be large below the 
barrier. In this case, the breakup yields are large but fu- 
sion is not much enhanced because the threshold anomaly 
is absent. This may explain why a subbarrier fusion en- 
hancement was not observed in ^^F-|-^°^Pb. Simplified 
coupled-channels calculations were performed to explore 
the effects of coupling to both inelastic excitations and 
transfer degrees of freedom on fusion. Furthermore, the 
radius of ^^F was adjusted in the calculations to study 
the change in the fusion excitation function. The re- 
sults suggest that fusion may be suppressed at energies 
below the barrier. In contrast, large subbarrier fusion 
enhancements were observed for fusion of the neutron 
skin nucleus ^He on ^"^Bi and ^"^^U. Further experiments 
are required to examine whether the differences observed 
between the ^He- and ^^F- induced fusion are due to 
breakup or other reaction mechanisms. Measurements 
using neutron halo nuclei, such as ^^Be and ^^Li, and 
proton halo nuclei, such as ®B and ^^P, would provide 
useful additional information. 
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